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ABSTRACT 

In  Part  I  the  Vlasov  equations  are  reduced  to  a  general 
form  for  solving  initial  value  problems.   These  equations  are 
then  studied  from  the  point  of  view  of  type  and  of  important 
scale  lengths. 

It  is  noted  that  transverse  flows  depending  on  only  one 

space  variable  depend  on  c/a3   for  small  mass  ratio  but  c/to  . 
—^ '^        '  pe  pi 

is  the  natural  length  for  all  other  cases.   Steady  flows  on  a 

scale  length  c/oo  .  are  studied  as  asymptotic  limits.   Estimates 

show  that  the  fields  and  currents  cannot  grow  faster  than  c/oo  A 

^  ^    pe 

where  A  is  Alfven  speed  unless  the  electrons  are  hot. 

In  Part  II  the  time-dependent  one-dimensional  equations 
with  transverse  field  are  investigated  by  various  methods.   The 
asymptotic  treatment  of  electrons  is  determined.   Correctly 
posed  problems  for  the  initial-boundary  value  problems  are  dis- 
cussed and  also  some  limiting  cases.   This  study  is  made  in 
Lagrangian  coordinates.   The  results  of  computations  for  weak 
Riemann  problems  where  there  is  no  breaking  and  rarefaction 
piston  problem  are  presented  along  with  the  methods  Involved  and 
justifications  for  certain  assumptions  such  as  neglecting  third 
moments. 

Finally,  a  model  for  strong  shocks  is  presented  in  Part  III. 
It  is  shown  that  under  rather  weak  assumptions  the  final  state  is 
uniquely  determined  by  the  initial  stage. 


Ill 


Introduction 

These  notes  represent  a  variety  of  results  all  based  on  a 
charge  neutral  collisionless  self -consistent  model  of  a  plasma 
of  two  species,  electrons  and  ions,  in  an  electromagnetic  field. 
These  requirements  mean  that  we  are  dealing  with  two  distribution 
functions  each  satisfying  Boltzmann's  equation  without  a  collision 
term.   The  forces  are  the  electromagnetic  forces  ruled  by 
Maxwell's  equations  assuming  the  speed  of  light  to  be  infinite. 
Thus  the  net  charge  is  zero  and  there  is  no  displacement  current. 

In  Part  I,  §1  the  equations  are  reformulated  to  lead  to  an 
algorithm  for  finding  a  solution  as  a  function  of  time.   The 
limiting  case  of  zero  density  is  briefly  considered  in  §3.   In  §4 
the  interesting  case  in  which  the  mass  ratio  of  the  two  particles 
is  regarded  as  asymptotically  small  is  studied  with  a  view  to 
determining  the  important  length  scales.   Vital  differences 
depending  on  the  geometry  of  the  fields  are  brought  out.   In  §5 
and  §6  asymptotic  properties  are  pursued.   The  most  interesting 
result  for  the  time  dependent  case,  §6,  is  that    in  the  zero 
pressure  case  we  must  have  density  proportional  to  field  magnitude 
on  each  field  line. 

In  Part  II  the  one-dimensional  problem  with  transverse 
field  is  considered  in  some  detail.   Actually  to  compute  this 
problem  one  must  use  asymptotic  theories  or  some  other  hypothesis 
to  make  it  tractable.   Various  possibilities  for  electrons  are 
considered.   For  ions  the  simplest  notion  is  to  discard  the  higher 
than  second  moments.   Some  computations  were  carried  out  and  in 


part  checked  for  their  error,  see  §2  for  the  general  results  and 
§6  for  the  detailed  comment  on  the  figures. 

In  §4  the  rather  delicate  question  of  what  is  a  correctly 
posed  problem  to  correspond  with  some  physical  problems  is  looked 
into . 

§5  and  §7  discuss  the  computational  method  for  the  differen- 
tial equations  and  the  third  moments  respectively. 

Part  III  hazards  a  possible  approach  to  the  structure  of 
strong  shocks.   This  may  very  well  involve  turbulence  which  is 
itself  so  complicated  that  we  cannot  study  strong  shocks  by  using 
it.   However  some  qualitative  properties  would  give  us  useful 
shapes  in  .principle  and  this  part  deals  with  a  shock  model  using 
a  minimum  of  special  properties.   It  is  a  model  in  two  parts.   The 
first,  §2,  is  a  quasi-structure  for  an  inner  shock.   The  second, 
§3,  is  for  the  exterior  overall  shock  structure  for  a  situation 
where  there  is  an  inner  shock. 

The  author  is  indebted  to  H.  Grad  not  only  for  helpful 
discussions,  but  for  offering  an  inspiring  scientific  point  of 
view  towards  this  subject. 


Part  I  -  General  Equations 

1.  The  General  Equations  of  Motion 

We  consider  a  general  two-component  charge  neutral  gas  of 
Ions  and  electrons  with  distribution  functions  f.(u,x,t),  charges 
e.  and  masses  m. ,  1  =  1,2,  respectively.   The  self -consistent 
equations  of  motion  are  then. 


^it  +"^'^^1  "^""l^i^^ +^^''^^^*^u^l  =  °  '     1  =  1,2  , 


curl  B  =  [i.J 


curl  E  =  B, 


q  =  0  ,    div  B  =  0 


where  E,  B  are  the  electric  and  magnetic  fields,  the  current  is 


J=Z/  e.f.u|du|  and  the  charge  q  =  2q. 


Ze.f.    du  . 
1  1 '   ' 


Introducing  the  length  /m^mp/]  e,  ep|  p*|j.  where  p*  is  some  ion 
density  and  normalizing  the  other  variables,  we  obtain  with  s-i  =  -So 


(1)    f .  ,  +U-V  f .  + 

^      '  it       XI 


-1 


(E  +uxB)-V  f .  =  0 

^         '   u  1 


with  1=1  use   e 


with  1=2  use  -£ 


where  e  =  Jm.-.  /mp  =  Jm   /m  and 


curl  B  =  J 


(2: 


■curl  E  =  B 


t 


^i  =  q2  '   ^-'-^  B  =  0 


with 


(3)  J  =  e  f   (fi  -f2)u|du| 


(1|)  q^  =  /  f.|du| 


From  the  equations  it  is  clear  which  quantities  are  vectors. 

With  these  equations  we  shall  have  certain  data  given  at 
t  =  0  which  we  do  not  yet  specify.   We  shall  also  keep  in  mind 
the  case  in  which  data  is  given  on  a  surface  '6  and  we  are  con- 
cerned with  solutions  in  its  exterior  or  interior. 

Our  first  object  is  to  reduce  these  equations  to  the  form, 

U^  =  A(U,V) 

(5) 

DV  =  G(U) 

where  A  is  a  non-linear  partial  differential  operator  in  the  space 
variables  and  an  integral  operator  in  the  velocity  variables;  D  is 
a  linear  differential  operator  on  V  with  coefficients  depending  on 
U  and  containing  no  time  derivatives.   D  also  is  invertible  pro- 
vided the  density  of  the  gas  does  not  tend  to  zero  anywhere. 

The  advantage  of  this  form  is  that  we  can  also  write  it  as 

t 
U(t)  =  U(0)  +  /  A(U,V)dt 
0 

V  -  D~-'-G(U) 

where  U(0)  are  specified. 


From  this  form  one  is  led  naturally  to  various  difference 
schemes  which  may  or  may  not  be  stable. 

This  method  was  developed  by  Morton  [3]  in  his  one- 
dimensional,  two-fluid  study  of  compression  waves  and  collision- 
less  shocks.   He  found  conditions  for  a  stable  difference  scheme 
provided  the  solution  existed  and  he  also  considered  cases  where 
only  weak  solutions  with  shocks  occur. 

We  note  that  our  equations  are  not  in  the  desired  form  (5) 
because  E  occurs  only  in  one  of  the  equations,  B,  =  curl  E,  and 
not  in  any  equations  without  a  time  derivative. 

The  desired  equations  are  found  by  taking  moments  of  (l) 
with  respect  to  the  velocity  variables  and  rearranging  the 
variables.   We  proceed  to  do  this: 

Integrating  (1)  with  respect  to  u  we  get  the  two  mass  con- 
servation equations: 


q.   +div  /  uf . I du|  =  0 


or  introducing  the  notation 

(6)  <Q>.  =  f  q{u)r  du  , 


(7)  J^  ^i  +  di^  <^>.  =  0 


Multiplying  (l)  by  u  and  integrating  we  obtain  the  six 
momentum  equations: 


(8) 


J^  <u>   +div  <uu>_  +1       _^l  (-Eq^-<u>^    X  B)    =    0    . 


Here  uu  denotes  the  tensor  u^u^,  k,.^  =  1,2,3. 

For  future  use  we  also  obtain  here  the  next  highest  moment 
equations  by  multiplying  by  uu.   Thus 

(9)   1^  <uu>.  +div  <uuu>.  +  I   _^|(-E  .<V^uu>  -<(u  xB).V^uu>^)  =  0  . 
Note  dlv<uuu>=  2  J^<u^u^u.>,  E.<V^uu>.  =  E^<up.  +  E/u^  >. 

u 

and  finally  <  (uxB)-V^uu>.  =  <(u  xB)^u^+  (uxB)^u^>.  . 

We   introduce   the   density   p   and   the   mean  fluid  velocity  U  by 


(10) 


(fl  +  e2f2)ldul    =<1>^   +<"^>2  ^    ^^^^^^  <^>i 


pU  =<u>^  +£^<u>2 


where    (2)    and    (4)    vjere   used. 

From   (2)    and    (3)    we   then  have   curl  B  =  < eu >   -    < eu >     or, 

introducing   the   mean   ion  and   electron  velocities 


<u>    =<1>U^   =    (l  +  e^)"^(pU+£   curl  B) 

TT          TT  ,        curl   B 
or     U-^   =  U  +  e 1 

(11) 

<u>     =  <1>U2   -    (1+  e^)~^(pU-  e~^    curl   B) 

,,        -1   curl  B 
or     U^   =  U  -  e 


From  (7),  (10),  (4)  we  see  that  there  is  really  only  one 
equation  for  conservation  of  mass: 


(12)  p  +dlv  (pU)  =  0  . 


From  (8)  and  (2)  we  obtain  by  linear  combinations  the  equa- 

2 

tlons  for  the  total  momentum  <(u)>  +  e  '\'u^  =  pU  and  the  current 

<'u^  -<Cuy     =    e~      curl  B: 

(13)  (P^U  +^i^  (<uu>   +£<uu>  -  curl  B  xB  =   0 
and 

(14)  ^    (curl   B)  +  e   div   (<uu>   -<uu>    -p(E  +  UxB) 
OL  12 

l-e^ 
+  curl  B  xB  =   0 

where  we  have  used  (11). 

After  some  calculation  we  find  from  (13)  and  (11) 


(15)       (pU)^  +  div    (pIKJ  +M^  +  E%^)  +  curl   B-V  -5^1^-1-  curl   B  xB   =   0 


where   M.    =  < (u  -  U.  ) (u  -  U. )  >  . 

Prom  (14)  and  the  equation  for  B,  in  (2)  we  obtai 


n 


B,  =  -curl  E  =  curl  (UxB)  -curl  —  (curl  B), 
t  ^  '  p  ^       '  t 

-  e  curl  p~   div  (<^uu^  -<^uu)>  ) 

1-e^      1 

-  curl  —  (curl  B  xB)  . 

e        p  ' 


Next  introducing  the  variables  P,  U"*  by 


pU*  =  curl  B 
(16) 

pF  =  B  +  curl  U* 

into  the  preceding  equation  we  obtain 

(17)*    (pF^t  "  ^"^^  [UxB  +  p~^(div  pU)U*  -  ep~^  div  (<uu>  -<uu>  ) 


Returning  again  to  equation  (l4)  and  using  (l6)  we  can  find 
an  equation  for  E  using  from  (2)  the  expression  for  B, .   This 
yields 

(18)*   E+  p"''"  curl  curl  E  =  -UxB+ep""^  div  (<uu>  -<uu>  ) 

+  Iz^  jj*  xB  . 


We  can  put  (17)   and  (l8)   in  better  form  by  using  along  with 
p,  +div  pU  =  0, 

div  (<uu>  -<uu>  )  =  div  (M-^^-  ¥\^)  -        g   pU*.VU* 


+  -  U-V(pU*)  +^  pU*.VU+£-  (div  U)U* 
Thus 


(17)   p(F^+U-VF)  =  B-VU-curl  (U-V)U*+  (U-V)  curl  U*  +  div  U  curl  U' 

.   2 

-  curl  (U*.VU) +^:-^^  curl  ((U*.V)U* 

-  (U*xB))  -  curl  -  div  (M^-Mg) 


and 

(18)  E+  p~^  curl  curl  E  =  -U  x  B  +-  div  (M-^  -  Mg)  -  ^^^^-   (U*.VU* 

-  U*xB)  +U-VU* +U*-VU+ p"-^  (dlv  pU)U*  . 
From  equation  (17)*  we  only  note  further  that 

(19)  1^  div  pF  =  0  . 

We  also  put  (15)  in  the  usual  form,  using  (12), 

(20)  l^+U-VU  +  p"-^  dlv  (M^  +  e^Mg)  +U*.VU*  -  U*xB  =  0  . 

We  note  a  few  interesting  facts  about  our  set  of  equations 
which  are  now,  (l)  for  f.  and  (12),  (l6),  (17),  (l8)  and  (20). 

I.  (1),  (12),  (17)  and  (20)  can  all  be  written  in  conservation 
form.   For  steady  one-dimensional  flows  mFrn  +  — dlv  (M-,  -Mp)  is 
constant  as  well  as  the  mass  flux  m  and  momentum  flux.   Here  sub- 
script T  indicates  transverse  component  with  appropriate  sign. 

II.  It  is  not  a  priori  clear  that  if  there  is  a  solution  of  (l), 
(12),  (l6),  (17)  and  (20)  that  it  is  a  solution  of  (l)  and  (2). 
This  can  be  proved  under  rather  general  conditions  for  one-dimen- 
sional flow. 

III.  The  highest  derivatives  of  U*  (essentially  the  electron  speed  ) 
that  occur  are  second  derivatives  in  (17)  and  these  occur  only  in 
the  curl  (U^.W)  term. 


The  equations  (1),  (12),  (l6),  (17),  (l8),  (20)  are  now  in 
the  desired  form  (5)  as  far  as  the  time  derivatives  are  concerned. 
The  variables  'U'  are  f.,  p,  pU,  F  and  the  variables  'V  are  B, 
\J*,    E. 

Initial  data:   From  the  general  form  it  is  reasonable  to  prescribe 
in  our  case  f . ,  p,  pU,  and  F  initially  with  p,  f.  compatible. 
However  F  is  a  somewhat  artificial  variable  and  we  can  prescribe 
instead  B  provided  div  B  =  0.   Then  from  (l6) 

1 
pF  =  B  +  curl  —  curl  B 

initially  and  hence  div  F  =  0  at  t  =  0.   Thus  from  (19) 

(21)  div  pF  =  0  . 

Physically  we  would  also  expect  to  require  curl  B  =  J 
initially.  We  shall  return  to  this  point  in  discussing  the 
computations. 

The  next  question  is  whether  we  can  solve  for  B,  U  and  E. 

Solving  for  B,  U*  and  E:   We  regard  F  and  p  as  given  in  (l6). 
By  (21)  we  may  write 

(22)  pF  =  curl  A 

and  A  may  be  regarded  as  known.   Let  ^  be  a  vector  and  ^   a  scalar 
that  together  satisfy 

-AT  =  p(A  -  T)  +  pV(|) 

(23) 

pA(t)  +  Vp  .  Vi|)  =  -div  pA  +  Vp.T 

10 


and  which  have  square  integrable  graldients  in  infinite  space. 

We  shall  show  that  the  solution  exists  for  p  j^  0  later. 
Meanwhile  we  set  B  =  curl  ^,  u*  =  A  -  Y+V(f),  and  find  that  (l6) 
is  satisfied  as  follows:   From  {2J>) , 

A  dlv  ^  =  p  div  ¥ 

and  hence  since  dlv  Y  — »■  0  at  oo  ,  div  ^  =  0  by  the  maximum  principle, 
Therefore  curl  B  =  curl  curl  Y  =  -A"?  =  Ap  -  p^  +  pVtj)  =  pU*  and 


curl  U*  =  curl  (A  -  ¥)  =  pF  - B 


as  required.   Furthermore  div  B  =  0  also.   Hence  the  problem  of 
solving  for  B  and  U*  is  reduced  to  solving  (23). 

If  (p,F)  — ^  (p  ,F  )  fast  enough  as  x  — »•  oo  in  one  direction, 
say  the  -x  direction  (or  upstream),  we  may  take 


A  =  xA,  +  yA^  +  zA,  +  A 


* 


A-,  ,  Ap,  A^  are  all  constant  vectors  and  A  — »■  0  rapidly  at  oo  .   We 
set 

Y  =  -(xA^  +  yA^  +  zA  )  +  ^*  . 

Then  we   find   that   our  equations    (23)    become 

(24)  AY*  -  pY*   =   pA*-  pV([) 

(25)  pA(t)  +  Vp  .  VcJ)   =   -dlv   pA*-dlv    (  pxA^+ pyAg+ pzA^) 

+  Vp.(-xA^  -yAg-  zA^)  +Vp.T* 

11 


or 

(26)  pA([)  +  Vp.Vcf)  =  g  +  Vp.Y* 

where  g  -*  0  as  fast  as  p  -  p^,  Vp  etc. 

These  equations,  (24)  and  (26),  can  be  solved  by  iteration 
if  p  is  almost  constant,  p  ~  p   ^  0  or  existence  can  be  established 

by  a  complete  continuity  argument  along  with  uniqueness  as  follows: 

2 

Let  us  first  prove  uniqueness  if  VB  lies  in  L   and  F  is  given. 

To  do  this  we  return  to  equation  (16)  with  F  =  0  and  find  that 

since  B-curl  U*  =  div  Q  -  )    \'^^\      ^^  ^^   integrate  over  space  we 

i      i 

find  /  ')    -  l^^l  dx  +  /  B-B  dx  =  0.   The  only  solution  of  the 
homogeneous  equation  is  B  =  0. 

Here  we  note  for  future  reference  that  Q  =  -div  B  x U*  and 
hence  for  the  uniqueness  of  boundary  value  problems  we  should 
require 

^  '  ■'    The  normal  component  B  x  U*  prescribed  on  surface  ^  . 

We  consider  next  integral  equations  corresponding  to  (24) 
and  (26)  restricting  ourselves  to  p  7^  0.   The  fundamental  solution 
for  the  operator  A  -  p   is  r~   exp  -  \/p~  r  and  for  A  is  r    with 
V  =    I  X  -  X  1  .   Thus 


^  =|1  (-  Vp^  ^^.  ^-  +  g)|dx|  =  L^(K^^)+f  , 
(28) 


W*      z= 


r'^   exp  (-/p^  r)((p  -  p^)Y*-  pV^  +  pA*)|dx|  =  L^i^,^*)  +  h 


12 


We  note  that  V(-)  and  r   exp  - /p^  r  both  lie  in  L^  at  co  . 
Thus  one  finds  with  \\^\\'^  =   j    |  V({)  |  ^  |  dx  |  and  ||T||^  =  /  (^*.Y*)|dx| 
that  by  classical  methods 

||VL^(<j),^*)||2  <   a(||(j)||2  +  ||^*||2) 

where  a  is  an  absolutely  constant  multiple  of  the  maximum  of  Vp/p 
and  P  -  Pq  and  P  is  bounded.   Hence  L-,  and  L^  combined  are  contrac- 
tion operators  for  Vcf),  Y  if  a  is  sufficiently  small,  that  is  the 
pair 

are  contraction  operators.   Thus  the  standard  iteration  scheme  will 
converge. 

Existence  also  follows  because  the  operators  are  completely 
continuous.   This  Is  easily  proved  by  standard  methods.   We  have 
already  established  uniqueness  so  now  we  have  existence  provided 
Vp/p  is  bounded. 

This  argument  also  works  for  some  boundary  value  problems. 

The  electric  field  E  can  be  found  by  Inverting  the  equation 
for  curl  curl  E  and  noting  that  dlv  pE  is  given  as  a  function  of 
the  other  variables. 


15 


2.  Reduction  of  the  Momentum  Transfer  Equations 


We  note  further  div  <uuu>=  ;§^  <^ ^ -i ^v^ »  / 


=  Jt  (^k.^  +  ^kU^  +  Vj^,+UAi''<VA>^  ^"^^^^  ^""^    subscript  i 

J 

has  been  suppressed  and 

Ki  =%  =<(^k-Uik)^^^-UiP>. 


where  U. .  are  the  three  components  of  U. .   We  then  expand  (9)  and 

recall  that  along  with  the  other  equations  we  must  obtain  an 

identity  if  all  the  second  and  third  moments  vanish.   Thus  we  obtain 

B  =  (B-,,B^,B^),  and  with  e,  .   vanishing  if  k  =  i,  i  =  r  or  k  =  r, 
^  1 '  2'    5  kir         °  ' 

and  equalling  ±1  as  kir  is  an  even  or  odd  permutation  of  1,  2,  J>, 


4r  ^:  n    +  ^-   (M.,  ,  +M.,U.,  +M.,  U.  ,  +M,  „U.  .) 
dt  ki        dx  .    jki,        j£    ik    jk  li    kl   ±j' 

-  U.  „   3^  M,  .  -  U.,  ^r^  M,  .  -  5(e.  .   M.  ,B^+  £,.   M.,  B^)  =  0 
i£  dx .   kj    ik  ax.      £^         ^    kir   i^  r    £±r      ik  r' 

where  M,  ^,  U.  etc.  is  for  ions  or  electrons  according  as  6  is  s 

1 
or  -  — . 
e 


3-  The  Limiting  Case  as  p  -»  0 

We  observe  that  if  p  -»■  0  at  a  finite  point  then  our  equations 
(24-26)  are  singular  since  V  log  p  =  Vp/p  becomes  infinite.   One 
needs  either  an  existence  theorem  or  to  establish  that  p  -*  0  only 
at  oo .   In  the  sheath  problem, the  simple  Rosenbluth  model  [l]  has 
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p  'v  yx  -  x^  and  this  coefficient  is  singular  but  in  Sestero's  models 
p  vanishes  only  at  oo ,  [2]. 

The  other  equations  are  all  non-singular  as  p  ->  0  and  in 
principle  give  no  difficulties,  except  that  the  Alfven  speed  tends 
to  00  .   We  note  only  that  U"^  must  be  bounded  so  curl  B  -*  0  but  that 
F  may  be  singular  with  pF  bounded.   We  see  this  by  considering  the 
equation  for  pF  instead  of  F. 

4.  The  Limiting  Case  of  e   Small  and  the 
Associated  Wave  Lengths 

Here  we  observe  that  there  are  two  distinct  cases  to  consider; 
We  set 

curl  Q  =  curl  (U*  x  pF  -  -  div  M^ ) 

and  the  two  cases  are 

Case  A:     curl  Q  =  0 
Case  B:     curl  Q.  y^  0    . 

This  crucial  expression  is  found  in  the  differential  equation 
for  F.   It  is  the  only  term  that  becomes  singular  as  e  -*  0.   As 
shown  in  (17),  this  term  is 


curl  Q  =  curl  ((U*-V)U*-  U*  x  B )  -  curl  |-  div  M2 


(29)  =  curl  (U""  X  (B-  curl  U  )  -  -  div  Mg) 

=  curl  (U*  X  pF  -  |-  div  M2)  • 
15 


Here  we  have  used  the  fact  that  Mp  is  0(— p)  if  the  pressures  of 

e 
ions  and  the  electrons  are  commensurable. 

We  first  make  an  estimate  for  how  large  the  currents  and 

magnetic  field  can  become: 

From  (17*)  using  (l6) 

B. (pF),  =  -  curl  ^  +  0(1) 
or  by  (16)  dotted  with  B, 


(4  IbI^),  -  B.curl  U*   =    -  -   B.curl  Q  +  0(1) 


On  integrating  over  space,  we  see  that  we  may  integrate  by  parts 
if  we  assume  that  all  gradients  vanish  at  large  distances.   Thus 


(50: 


^J     '2 


^  1b|^+  I  curl  B|^]|dx|  -U*.div  M2 1  dx  |  =  0(1) 


since 


B'curl  Q  dx 


curl  B- (U*  X  F  -  -  div  Mo)|dx| 

p      2  '   ' 


pU*. (U*  X  pF  -  -  div  Mp) Idxj 


=  £    U*.div  Mo  dx 


and 


B' curl  U, I dx 


r*   TT* 


pU  .U^ldxl  =  -  ^ 


^   ^p|U*|^|dx|  +0(1)  . 


The  term  involving  the  electron  pressure  tensor  represents  of 
course  the  work  done  on  the  electrons. 
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Equation  (JO)  asserts  that  provided  U*.dlv  Mp  =  o(— ),  i.e. 
the  work  done  by  the  stresses  in  the  direction  of  the  current  is 
0(1)  the  fields  remain  finite. 

On  the  other  hand  we  see  that  E  ~  0(— )  unless  Q  itself 
vanishes.   The  quantity  pF  is  a  modified  field.   If  the  gradients 
of  B  are  small  pF  'u  B.   Hence  if  Q  vanishes  we  have  a  particular 
configuration  somewhat  like  a  force-free  field. 

Thus  it  is  to  be  expected  that  there  will  be  two  kinds  of 
flows  —  those  in  which  the  term  in  curl  Q  dominates  and  those  where 
it  is  small. 

The  natural  wave  lengths:   The  two  natural  wave  lengths  are  for: 

Case  (A):   The  basic  unit  used  here  for  normalization,  i.e., 
c/cD   .   If  we  let  e  — ••  0  and  drop  e  compared  to  1  we  have 

p,  +  div  pU  =  0 
U^  +U-VU  +  p'^  div  (M^  +  e\^)^  u^'-VU*-  U*  x  B  =  0 
p(F,  +U-VF)  -  +  B-VU  -  curl  (U-V)U*+U-V  curl  U* 
(31a)    ■  +  (div  U)  curl  U*-  curl  U*.VU  +  curl  |-  div  M^ 


P 


E  +p"-^  curl  curl  E  =  -U  x  B  -  |-  div  M2  +  |  U-VU* 

+  -   U*-VU+  p"^  div  (pU)U* 


pF  =  B  -  curl  U*  ,    pU*  =  curl  B  . 
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A  special  case  which  has  good  asymptotic  properties  occurs 

if  curl  (U*.VU*)  ^   0  but  curl  (U""  x  B+e^  div  ¥ir^)    =   0.   We  intro- 

1/3 
duce  as  a  new  length  scale  e  ^^x  and  in  this  unit  obtain 


p,  +  div  pU  =  0 
U^  +U-VU  +p~-'-  div  (M^  +  e^Mg)  -  ^^^^   ^  x  B  =  0 
p(F,  +U-VF)  =+B-VU+curl  (.^HlL^  .  v  ^H^i_^)  _  curl  |-  div  (M.  -  Mp) 

E  +U  X  B  =  ^  div  (M^  -  Mg)  -  £H^.  V  ^.^  -  B-2/5  (£H£LB  ^  g) 

and 

pF  =  B  . 


The  advantage  of  this  model  is  that  the  electrons  are  truly  adia- 
batic  and  the  ions  see  mainly  electric  field. 

Case  (B):   The  basic  length  (and  time-scale)  should  balance 
the  two  sides  of  (17).   Hence  t  and  x  should  be  scaled  up  by  a 
factor  e,    i.e.,  x  =  ex  ,  to  the  ion  length,  c/co  .  .   The  equations 
reduce,  on  dropping  the  primes  again,  and  neglecting  terms  of  order 
e  to 
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•p.    +  dlv   pU  =   0 


U^  +  U-VU+p    -^    div    (M-^+  E^n^)    -    ^^^^   ^    X  B   =    0 


(31b)  pU*   =    e   curl   B    ,         pF   =  B 

2 

E  +U   X  B   =  ^^  div  M^    -    ^^^-^   ^   X  B 
P  2  p 

2  -,    „ 

B.    +div   UB   -  +B-VU+  curl    (—  div  M^    -    ^^^     ^   x  B)    . 

t  P  £^  p 

Note  that  the  electric  field  is  of  order  l/e  in  Case  (A) 
unless  Q  =  0  but  is  of  order  1  in  Case  (B). 


5.  Motion  of  Particles 

Case  (A):   The  ions  are  governed  by  the  electric  field.   The 
electrons  may  be  treated  adiabatically  only  for  thermal  speeds  less 
than  order  l/e.   This  means  that  the  corresponding  P-factor  for  the 
electrons  must  be  small.   There  is  a  special  adiabatic  invariant. 
Strictly  speaking  Q  ~  U  x  pF. 

This  case  occurs  for  example  in  studying  the  structure  of 
low  p>  transverse  shocks,  where  there  is  only  one  space  variable. 
See  Morton  [5],  Auer,  Hurwitz,  and  Kilb  [k],    Kilb  [5],  Morawetz 
[6],  Marder  [7].   See  also  Gardner,  Goertzel,  Grad,  etc.  [8], 
Morikawa  and  Gardner  [9].   In  fact,  this  is  the  appropriate  length 
scale  for  any  field  perpendicular  to  the  axis  of  a  single  inde- 
pendent  (space)  variable.    In  general  for  two  space  variables  even 
in  a  transverse  field  oo  .  is  the  scale  length.  Case  (B). 


This  length  first  appeared  in  the  Adlam-Allen  wave  [10] 
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Case  (B):   Here  E  =  0(l)  in  e  and 

2 

E-B  =  —  B-  dlv  (M^-  Mg) 

which  as  we  shall  see  may  be  0(1).   However  we  may  treat  the 
electrons  adiabatically  even  if  their  pressure  is  finite  for  the 
following  reasons  (despite  the  usual  requirement  E-B  =  0(e)): 

Suppose  the  pressure  is  finite  since  the  mass  is  small.   Most 
electrons  must  have  large  velocities,  i.e.,  to  have  e   /  fu  |du| 
finite  and  /  f|dul  also  of  order  1  we  must  have  u  =  U/e  where  U  is 
finite. 

For  the  length  scale  of  Case  B,  c/co  .  ,  the  equation  (l)  for 
the  electrons  becomes 


e^f^  +£U-Vf  -  (eE  +  U  x  B)-Vyf  =  0  . 


This  has  the  same  asymptotic  behavior  as 


U-Vf  -  ^  (eE  +  U  X  B)-Vyf  =  0  ; 


the  equation  describing  particles  in  a  magnetic  field  B,  electric 

field  eE  and  a  mass  to  charge  ratio  e.   Since  eE-B  =  0(e)  we  can 

2 

use  the  standard  adiabatic  theory  [15]  and  thus  e  M^  is 

/  P   0   0 
0    p^  0 
\0    0    p^ 

in  the  framework  of  the  field  line  as  x-axis  with  p^^  =  p^/|b|  . 
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p^  Is  constant  on  the  field  lines.   To  order  e  the  energy  is 

conserved  for  each  particle  and  thus  p  =  3p  (1 ^— ?^) .   Note 

that  since  eTvi  =  0(£  )  we  have 

E-B  ~  -  B-div  e^Mo 
P  2 

•V  —  B  •   -rg-  ^r—   R    1  -   ^ 

p     R  ^   (^    ^|g|2 

where  R  is  the  radius  of  curvature  of  the  field  line  and  s  is  arc 
length  on  it. 

This  situation  occurs  for  shocks  in  oblique  fields,  see 
Morton  [5],  Drummond  [11],  and  also  Kever  and  Morikawa  [12].   It 
also  occurs  for  a  transverse  shock  with  finite  p.   The  one- 
dimensional  solution  to  this  problem  requires  a  discontinuity  to 
admit  a  shock.   R.  Morse  [13]  has  examined  the  ion  stability 
numerically  in  a  related  model  and  found  evidence  of  cellular 
behavior  ahead  of  a  shock. 

One  should  also  mention  the  original  wave  of  this  type,  the 
Saffman  wave  [l4] . 

6.  Steady-state  Solutions  near  Limit  e  -»  0 
Case  (B).   Ion  Length 

We  investigate  some  asymptotic  properties  of  the  limiting 
case  in  the  steady  state.   The  underlying  mathematical  question 
is:   Wliat  are  the  asymptotic  connections  between  our  original 
equations  and  the  limiting  case  (B)  described  by  ( 31 ) .   We  pose 
this  only  for  an  equilibrium  state.   We  are  looking  for  the 
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properties  of  solutions  whose  natural  length  scale  is  the  ion 
length  c/o)  .  . 

P  -L 

Preserving  the  correction  terms  in  the  many-dimensional  case 
we  have  for  an  equilibrium  solution  from  (20)  with  U  =  0: 

p~^  div  (M^+  e^Mg)  +U*.VU*  +  U*  x  B  =  0 
pU*  =  curl  B  ,    div  B  =  0  . 

Since  U  -VU*  is  of  higher  order  we  see  that  the  characteris- 
tics are  no   longer  the  field  lines  as  they  are  if  this  term  is 
missing.   Instead  the  characteristics  satisfy  U^.dx  =  0  or 
curl  B'dx  =  0  and  hence  are  in  the  direction  of  the  current. 

We  examine  the  case  of  zero  pressure 

U*.VU*  +  U*  X  B  =  0 
div  B  =  0  . 


Scaling  up  to  the  ion  length  c/co  .  introduces  the  factor  e 
and  we  have 

U*.VU*  +  -^  (U*  X  B)  =  0 

£ 

pU*  =  curl  B  ,    div  B  =  0 


where  U*  has  also  acquired  the  e-factor. 

If  we  regard  U*  as  a  velocity  the  first  equation  is  for  the 

motion  of  a  particle  in  a  steady  magnetic  field  and  its  energy 

2  2 

U*'  is  a  constant  U*  along  each  current  line.   There  is  an 
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adlabatlc  Invariant  If  e  is  small  which  is  U*  -  (U*.B/|b|)  /|b|. 
U  consists  of  a  component  parallel  to  the  magnetic  field  which 
is  varying  with  the  field  derivatives  and  a  component  perpendicu- 
lar in  a  direction  that  moves  with  period  0(— )  around  the 
magnetic  field  line.   That  is,  if  the  second  derivatives  of  B 
are  bounded  as  e  — »•  0. 

(32)      U"  =  (Uf  -  |B|  ^x)^/^  B   ^  (|3|  ^)V2  .  ^^^^^ 


where  j  =  sin  —  k  +  cos  —  ^;  k,  i   being  a  pair  of  unit  orthogonal 

vectors  perpendicular  to  B,  ^   =        JBJds  along  the  field  line; 

U*,  u  are  constants  on  each  field  line, 
o' 

The  Limiting  Solutions  as  £  -»  0  with  Finite  Second  Derivatives 

Assume  there  is  a  limit  with  B  having  finite  second  deriva- 
tives as  e  -^  0.   The  adiabatic  form,  (52),  may  be  regarded  as  a 
differential  equation  for  B  of  the  form 

curl  B  =  pg(B)  . 

If  there  is  a  solution  to  this  equation  then,  since  g  con- 
tains the  rapidly  oscillating  vector  j,  some  of  the  second  deriva- 
tives of  B  will  tend  to  infinity  as  e  -^  0  unless  the  amplitude 
factor  of  the  unit  vector  j  tends  to  zero.   Thus  in  (32),  la  ~  0 
for  e  small.   Hence  U*  ~  U*B/1b|  with  U*  constant  on  each  field 
line . 

On  the  other  hand  div  pU*  =  0  or  B-Vp/|B|  =  0.   Hence  p/ | B | 
is  constant  on  each  field  line  in  the  limit  as  the  Larmor  radius 


0  if  the  pressure  is  zero. 
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Note  that  if  we  proceed  to  the  asymptotic  limit  directly  we 
obtain  no  condition  on  p  only 


or 


U*  X  B  =  0 


U*  II  B 


subject  to  div  B  =  0.   Hence  there  are  many  solutions  of  this 
limit  which  are  not  true  limits. 

It  seems  unlikely  that  this  solution  will  be  stable  since  it 
Is  embedded  in  a  family  of  solutions  (32)  (with  \i  not  necessarily 
zero)  which  are  highly  oscillatory. 

Let  us  consider  next  the  case  of  finite  pressure.  Intro- 
ducing  a  fictitious  electric  field  E  we  see  that  the  new  equi- 
librium equations  are 

E  -  -  -  div  (M^  +  e'^m^) 


U*.VU*  =  \   (E+  U*  X  B) 


div  B  =  0  . 

In  order  that  U  =  0(1)  one  can  show  that  E-B  =  0(e).   Hence 
the  pressure  tensor  is  restricted  to  satisfy 


B-div  (M^  +  e^Mg)  =  0(e)  . 


Then  with  the  guiding  center  interpretation  we  find  that  in 
order  to  have  a  solution  of  0(1)  in  e, 
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„*      E  X  B  ^,   B 


B 


B 


,  -,  2  TT*2 

where   A   =  U 
o 


2 

E  X  B 


B 


where  U   is  constant  along  the  current  lines, 
o 

The  condition  div  curl  B  =  div  pU*  then  leads  to 


or 


B 


.V  (M_)  _  div  £^^ 


B 


B-V  (^^)  -  curl  (pE) 


B 


B 


+  curl 


B 


(-^)  .  ^^  =  0 


B 


B 


B 


Thus  we  have  a  differential  equation  for  the  pressure  tensor 
pE  =  M-,  +  e  Mp.   Or  we  can  regard  this  equation  as  a  restriction 
on  p  corresponding  to  the  zero  pressure  case. 
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Part  II  -  One-dimensional  Motion  with  Transverse  Field 

1 .  Equations 

We  consider  in  greater  detail  and  for  numerical  study  the 
case  in  which  the  only  independent  variables  are:   x  and  t  and 
the  field  B  is  in  the  z-direction.   The  dependence  of  f.  on  the 
third  velocity  component  drops  out. 

The  set  of  differential  equations  reduces  with  U  =  (u, v) , 
by  noting  that  curl  B  has  a  component  only  in  the  y-direction,  to 


Pt  ^  ^P^^x  "  ° 


u,  +  uu  +  p  -^  ^  (M,  +  e^Mo)       +p""^B^B^   =  0 
t    X   '^    dx  ^  1      2'x  comp   ^    3  3x 


V,  +  uv^  +  p  -^  ^  (M-,  +  e^Mo)       =  0 
t    X      ^        dx  ^  1      2^y  comp 


(I)    F^,  +uF,  =  -  ep   ^  -  :e-  (M.  -M^) 
'     3t     J)X  ^        dx  p  dx  ^  1    2^ 


1  a  1  a 

p  ^  '"  i  ' '-^'y  comp 


-U2X  =  (^  Bjx^x  -   ^3  -  P^y    U^  =  U;  =  F^  =  F^  =  0,  B2  =  B^  =  0 

2 

E  +  U  XB-  I  E^^^  =  ep"-^(M^  -  M^)^  +  uU%  p"^  (pu)^U*+  ^^  (U*xB) 


and  the  original  equations  (l)  for  f.  reduce  to 


(II)     f.,+uf.   +(E+uxB)-Vf  =  0,      i  =  l,2 
it     IX  "^   u  ' 


We  note  that  U-,  =  XJg   in  the  x-direction. 

Thus  the  momentum  equations  (32)  for  m|-,  =  N+,  M^p  =  M^, 
Mpp  =  L^  reduce  to 
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N±t  +N+xU  +3N^u^  -  5^2BM^  +  N^^^   =   0 
(III)         M^^  +M^^u  +2M^u^+N^V^-  5^{BL^-pFN^)   +1^^-^^^   =   0 
^±t  +  ^^±^)x+  2M+v^  +  6^2pFM^  +N^2±x   =   ° 


where  the  third  moments  are 


and 


N^+  =  /  f^(u  -u)^  dudv 


Ng-j^j.  =  /  f^(u-u)   (v-v^)  dudv 


^12+  "  /  f±(u-")  (v-v_j.)^  dudv 


6^  -  e  ,    5_  =  -  e  ^ 


Finally  we  will  want  to  know  the  electric  fields  which  are 
found  by  (l8 )  to  be 

2 

E^  =  -BV+  e  p-^(N^-  N_)^  -  ^^  P"^(B^)^ 

(IV) 

^2x  =  -  ^t  • 

Theoretically,  we  could  now  set  up  a  numerical  scheme  for 
solving  these  equations  in  the  following  way.   We  go  forward  in 
time  by  using  a  difference  scheme  on  (I)  and  (ill)  to  obtain  new 
values  of  p,  u,  v",  F,  L_^,  M^,  N_^  using  the  old  values  of  p,  u,  v", 
F,  L^,  M^,  N,,  B,  and  the  third  moments  differenced  in  x   to  obtain 
space  derivatives.   The  new  values  of  B  are  then  obtained  by 
solving  a  second  order  ordinary  differential  equation  for  F  (the 


27 


last  equation  of  (I)).   New  values  of  the  third  moments  are 
obtained  by  following  the  characteristics  of  the  two  transport 
equations  for  f,  i.e.  following  the  particles  and  estimating  the 
moments.   To  do  this  one  needs  the  electric  fields  which  are 
obtained  from  (IV).   However,  this  scheme  is  unrealistic, 
principally  for  two  reasons.   Firstly,  when  e  has  its  physically 
significant  value  approximately  I/60,  one  sees  from  (III)  that  we 
would  need  time  steps  that  are  small  compared  to  l/e  but  we  know 
that  significant  waves  will  take  place  on  the  given  time  scale 
of  order  1. 

The  simplest  way  to  resolve  this  dilemma  is  to  set  the  third 
moments  of  the  electrons  equal  to  zero.   This  can  be  completely 
justified  if  the  adiabatic  approximation  is  valid,  that  is,  if 
the  electron  pressure  is  small,  i.e.,  if  the  thermal  velocity  is 
0(1)  or  even  <<   l/e.   In  this  case,  see  [6],  the  distribution  is 
symmetric  in  u  -  u,  v  -  v,  to  lowest  order  in  e  and  not  only  may 
N^  ,  Np^   and  N, p  be  neglected  but  also  M_  is  small. 

From  (II)  for  the  electrons  we  see  that  for  e  small  we  must 
have  to  lower  order 

(1)  BL  -  FpN  -  0  . 

Using  (1)  and  eliminating  M  from  the  first  and  third  electron 
moment  equations  we  obtain  a  differential  equation  for  L   (or  N_ ) . 
We  find  after  using  the  conservation  of  mass,  the  first  equation 
of  (I), 

(2'  (|^-"Ie)^  =  ° 
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5  2-1  dx    — 

or  BV^L  F    is  constant  on  a  fluid  path  ^  =  u.   Using  (2) 

in  the  first  electron  moment  equation  we  obtain 
(3)  eM_  =|!(|^^-II^,(p-1l_)  . 


Note  that  we  need  £M_  to  find  P  in  (Ic). 

In  what  follows  we  shall  assume  that  in  any  event  we  can 
neglect  the  third  moments  of  the  electrons,  see  Appendix  A.   We 
shall  either  use,  from  (l),  (2),  (j),  assuming  e  is  small 

(4)        (Bp"5f-^)^/^L_=  Nq  ,    (FB~^p"^)^/V  =Nq 


where  N  is  constant  on  dx/dt  =  u,  and 


o 


(5)  BM_  =  ■|f""\(^+u^)(FB-V^)1/2 


or  we  shall  use  the  differential  equations; 


N    ,   +  N  ^u  +  3N  u^  +  2e    ^BM      =   0 
(6)  M.+M      U  +  2MU+NV+  e"-'-(BL    -  pFN    )    =   0 

L    ,    +(L  u)     +  2M  V     -  2e~-^pFM      =   0    . 


The  second  problem  is  connected  with  finding  the  third 
moments  for  the  ions.   It  is  reasonable  to  consider  only  problems 
in  which  the  third  moments  vanish  initially.   A  power  series 

expansion  in  time  shows  that  the  third  moments  develop  very  slowly 

u 
—  in  the  limit  e  — ►  0  like  t  .   On  the  other  hand,  it  is  very 

difficult  to  compute  their  derivatives  accurately  by  following 
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particles.   We  have  therefore  set  the  third  moments  equal  to 
zero  in  solving  (l)  and  (ill)  but  as  a  check  computed  their 
values  from  the  resulting  fields. 

The  equations  have  been  solved  with  the  third  moments  of  the 
ions  set  equal  to  zero.   See  Chew,  Goldberger  and  Low  [l6]  for  a 
similar  model. 

In  all  cases  where  the  differential  equations  with  the  third 
moments  set  equal  to  zero  appear  to  have  a  continuous  solution 
the  third  moments  are  so  small  that  they  cannot  have  any  qualita- 
tive effect.   In  cases  where  the  solution  is  about  to  break  (as 
the  zero-pressure  solitary  wave  breaks  if  the  particle  velocity 
goes  to  zero)  it  appears  that  the  third  moments  acquire  an 
oscillation  but  the  magnitudes  of  the  moments  will  not  alter  the 
situation  and  stop  the  breaking. 

2.  Summary  of  Numerical  Observations 

From  (I)  we  see  as  expected  that  the  'charge  separation' 
electric  field  E-,  is  0(l/e).   Hence  the  motion  of  the  ions  is 
dominated  by  this  field.   If,  in  (II),  we  set  the  third  moments 
of  the  ions  equal  to  zero  we  find  since  M_j_  and  (BL  -  FpN  )  are 


small 


{^+P~^)^+^i^+p'^)y.   =   0 


where  we  have  used  the  conservation  of  mass  p,  +  (pu)   =  0  from 
(Ic). 
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In  fact,  this  checks  with  our  numerical  calculations 
remarkably  well.   A  similar  argument  which  yields  a  different 
pressure  tensor  L_|_,  M_^,  N  can  be  carried  out  for  an  axial  pinch, 
see  Marder  [7].   His  results  also  involve  different  moments  for 
the  electrons. 

It  is  clear  that  dropping  the  third  moments  in  (II )  should 
be  a  good  approximation  if  the  thermal  speed  of  the  ions  is  small 
compared  to  the  other  velocities  in  the  problem.   What  is  remark- 
able is  that  even  with  thermal  speeds  of  the  order  of  .5  which 
corresponds  to  p  =  .5a  rarefaction  wave  which  cuts  B  in  half 
seems  to  satisfy  these  conditions  very  well,  see  Figures  yj , 
42,  43.   Then  the  largest  third  moment  at  time  T  =  4  is  still 
only  .02.   Smaller  values  of  p  are  still  better. 

On  the  other  hand,  if  we  generate  a  compression  wave  that 
does  not  break  in  the  case  where  we  set  the  third  moments  equal 
to  zero  we  find  that  the  possible  strength  of  the  maximum  com- 
pression declines  very  rapidly  with  pressure.   For  example,  if 
the  compression  is  induced  by  a  five  percent  change  in  magnetic 
field  the  thermal  speed  must  be  less  than  l/4  of  the  Alfven  speed, 
If  there  is  no  initial  magnetic  field  gradient  but  a  pressure 
gradient  (no  density  gradient)  which  induces  compression,  the 
wave  breaks  with  thermal  speeds  of  about  .2  or  P  =  .08.   For 

example,  the  neglected  third  moments  in  case  initially  p  =  . l4 

_4 
are  approximately  10   at  their  maximum,  i.e.  an  order  of  magni- 

h, 
tude  smaller  than  the  u  ,  where  u^^  is  thermal  speed. 

Let  us  recollect  here  that  when  a  wave  is  'saved'  from 

breaking  by  viscosity  two  properties  help:   1)  the  amount  of 
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viscosity  and,  li)  the  fact  that  this  coefficient  multiplies  a 
higher  derivative  and  hence  the  restoring  effect  is  amplified  as 
the  wave  tends  to  steepen.  Here  the  third  moments  occur  with  the 
same  order  of  derivative  as  the  second  and  the  magnitude  does  not 
grow  adequately.  Furthermore  any  apparent  oscillation  (it  is  not 
clear  that  this  represents  some  real  fluid  behavior)  will  have  an 
alternating  effect  on  the  breaking. 

3.  Equations  in  Lagrangian  Coordinates 
with  Neglected  Third  Moments 

We  introduce  the  Lagrangian  variables  T,  |by-^  =  -^  +  u^, 
-^  =   —  •^— .  The  differential  equations  simplify,  with  V  =  l/p,  to 

5V  _  5u 
'St      M 

2 

(7)  11=  -ifjtM^.Mj 


^  =  BV  -  F  . 


The  moment  equations  are  with  6=1,  5_  =  -e   , 
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(8,9)         ^V%  =  5^(BV%-FVNJ ±2  11 

1+e 

N  2M,   p.- 

1+e 

For  the  electrons  alternatively  for  low  temperatures,  (^,5), 

(10)  L_  =  g"/"(e)Fl/W5/2  ,    N_  =  ^  =  el/2(|)V-5/=^Bl/2F-l/2 

2 

eM_  =  |p  (VL_)^  . 

The  electric  fields  satisfy 

2  f^F 

(11)  E^  =  -BV  +  ^  (e(N^-N_)  -  (i^)B2)  ,    ^._v||  +  u||. 

^.  Correctly  Posed  Initial  Boundary  Problems 
and  Numerical  Method 

Simply  by  counting  we  observe  that  for  the  full  set  of  equa- 
tions (I)-(III)  we  may  prescribe  f+  initially  along  with  B.   On  a 
boundary  x  =  h(t)  we  may  prescribe,  if  the  region  of  interest 

corresponds  to  x-h(t)  >   0,  the  value  of  B  and  the  values  of  f^ 

dh 
for  u  >  ^.   If  the  region  corresponds  to  x  -•■  oo  we  require  B  to  be 

CI  w 

bounded  at  co  and  we  may  prescribe  f^   for  u  <  0,  x  =  oo  . 

In  making  approximations  by  discarding  the  third  moments  and 
using  the  differential  equations  for  all  second  ones  we  can  pre- 
scribe initially  the  six  moments,  V,  u,  v,  and  B  and  on  a  boundary 
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I  =  constant  we  can  prescribe  B  (or  one  other  quantity).   The 
same  is  true  if  we  use  the  asymptotic  system  for  the  electrons, 
We  have  reasoned  on  the  "basis  of  counting  alone  but  the  charac- 
teristic speeds  for  the  full  system  of  the  eleven  equations 
indicate  that  more  conditions  must  be  imposed- 

To  obtain  the  characteristic  speeds  from  this  cumbersome 
system  we  rewrite  the  equation  for  F  as 


and 


Fiji  -  Gj 


aG^  =  G  -  (M^  -M_), 


and  consider  the  limiting  case  as  a  — *■  0.   This  way  we  may  miss 
speeds  but  on  the  other  hand  we  do  not  appear  to  Introduce  any 
extraneous  ones.   Since  roots  are  multiple  anyway  we  cannot  use 
standard  characteristic  theory.   But  we  can  get  some  feel  of  what 


is  required.   The  characteristic  speeds  are  then  0,  oo  ,  ±/N  +  e  N_. 

This  suggests  that  we  ought  to  prescribe  some  conditions  on 
the  boundary  and  some  recent  work  indicates  that  the  appropriate 
quantities  are  M+,  B  and  u.   The  last  one  falls  out  if  we  look  at 
the  fluid  dynamic  limit,  B  =  0  and  the  condition  on  B  is  seen  from 
its  equation,  see  also  Marder  [7l« 


An  example  of  this  is  given  by  the  system  u,  =  Pxx'  ^t  ~  ^' 
Rewriting  this  as  a  limiting  case  of  u,  =  q  ,  aq  =  p  -  q,  p,  =  u 
with  a  -»•  0  we  find  the  three  speeds  0,  0,  0  from  the  formal  matrix. 
But  written  as  u,  =  q  ,  q,  =  u  ,  p  =  u  we  obtain  ±1,  0,  as  speeds. 


1>^ 


We  consider  next  correctly  posed  problems  for  the  low 
temperature  case.  I.e.  where  we  use  (4,5)  for  the  electron  moments. 

If  we  use  low  temperature  equations  then  the  equation  for  F 
may  be  written  as 

To  solve  this  equation  along  with  the  others  we  might  con- 

2 
sider  the  last  term  as  small,  i.e.   e  L   is  small  and  continue 

using  our  singular  perturbation  as  we  did  in  deriving  the 

expressions  (4,5).   Thus  we  set  F  =  F  +  F,  ,  B  =  B  -HB-, 


(12) 


o + 


J 


S^B  h\ 

pr^  =  BV-F   ,    ^  =  B,V-F,  . 

hi^  o  o    '         ^^2      1     1 

These  equations  replace  the  equations  for  F  and  B  in  (l) 
Combining  them  again  we  get 


(15) 


a^p 

dT      ■ 

^e 

5f 

0 

S^M,^ 

6T 

h' 

8^B 

BV-  F 

S^B. 

0 

=  B  V  -F 

O               0 

2F 


o   B  V^  T  J^^ 
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The  computations  for  low  pressure  electrons  were  carried 
out  using  a  variation  of  these  equations  which  will  be  discussed 
later.   See  also  Marder  [7]  for  a  better  scheme. 

Another  possibility  is  to  examine  the  system  directly.   We 
set 


2F   '2F    2B    2V 

2 

^p-  -^   log  (s  VL_) 


then  by  (7) 


F^  =  P.g  -  eM^.g    and    P..  -  eM^gg  =  F  (^  +  -^  +  -S^^— )  . 

Next  set 

(1^)  B^  =  R 

so  that  from  (7)  for  B  differentiated,  (7)  for  F  and  (7)  for  V  we 
have 

(15)  R^^  =  RV+BU^+P^^  -^^+ii 
and  also 

(16)  P^^„P(R.2!l  +  ^,.,„^^^  . 

We  now  compute  the  characteristic  speeds  as  limits  for  these 
equations  with  the  variables  R,  P  in  addition  to  V,  U,  E,  F,  v,  B, 
N+,  M^,  L^. 
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The  possible  speeds  are  ±Ji^_^+   e^N_ ) ,  0  (six  times)  and  oo 
(twice).   Thus  two  data  should  be  prescribed  on  ^  =  constant. 
One  of  these  could  be  u.   For  the  other  we  might  choose  M_  on 
^  =  constant;  say  M_  =  0  which  Implies  that  a  certain  symmetry 
prevails  at  the  edge.   This  would  be  resolved  in  a  boundary  layer 

of  scale  length  £   determined  by  the  coefficient  of  P  in  (l6), 

2      1  /2 
£   =   E    (VL_)  ^  /2F ,    i.e.  by  the  thermal  speed  of  the  electrons  with 

the  factor  e.   This  is  essentially  the  thermal  speed  of  the  ions 

if  the  two  pressures  are  both  roughly  the  same.   Ahead  of  the 

boundary  layer  one  may  reasonably  expect  the  asymptotic  behavior 

described  by  (13).   Recall  that  it  is  consistent  with  low  pressure 

electrons  to  neglect  the  third  moments  and  the  boundary  layer  both 

at  I  =  constant  and  at  t  =  0. 

Breaking  Waves 

The  computations  show  that  if  we  assume  our  adlabatlc 
formulas  for  the  electrons  or  neglect  the  third  moments  only  weak 
compression  waves  do  not  break.   (Rarefaction  waves  yield  small 
third  moments  and  are  quite  consistent.)' 

5.  The  Computational  Method 

A  grid  of  points  (l,J)  1  =  0...N,  j  =  0. . .M  is  introduced 
in  i,T  space  and  the  partial  differential  equations  (7),  (8)  and 
(9)  If  used  are  solved  in  the  following  way.   Every  T  derivative 


All  the  computations  described  in  this  and  subsequent  sections 
were  programmed  and  carried  out  by  Antoinette  Wolfe.   The  machine 
was  the  CDC  6600. 
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is  replaced  by  the  difference  D  which  satisfies 

2kD  U(l,j)  =  U(i,j+1)  -U(l,j-1)  and  the  first  order  ^  derivative 

by  d  where 

2hd  U(i,j)  =  U(i+l,j)  -U(i-l,j)  . 

Except  in  the  equation  for  B  the  second  order  derivative  is 
replaced  by  d*  with 

4h^d*  U(i,j)  =  U(1+1,J) +U(i-l,j) -  2U(l,j)  . 

At  the  edge  ^  =  constant  the  values  of  M  -  M_  needed  in  (7) 
are  extended  by  interpolation.   Thus  we  can  proceed  one  time  step 
forward  in  a  regular  way,  i.e.  all  variables  except  B  evaluated  at 
j+1,  1  =  0...h  can  be  computed  from  their  values  at  (i,j),  (i,j-l) 
for  all  1.   To  solve  for  B  we  use  the  method  of  Rlchtmyer  [3]-      The 
only  difference  is  that  P  =  1  in  his  problem.   The  method  makes  use 
of  the  fact  that  the  solutions  may  be  written  as  C(B,  -1)  +D  where 
B,  -1  satisfies  the  corresponding  homogeneous  equation  and  tends  to 
1  at  infinity  and  D  satisfies  the  Inhomogeneous  equation.   One 
prescribes  B  at  the  left  hand  mesh  point  to  determine  C. 

To  solve  the  adiabatic  equations  (10)  instead  of  (9)  we  use 
a  modification  of  (13).   At  each  time  step  we  solve  for  P  assuming 
M  =  0.   We  then  use  this  new  value  of  P  to  compute  M  from  (10) 
replace  in  the  equation  (7)  for  P  and  obtain  a  corrected  value  of  P. 
This  is  probably  unstable  if  the  pressure  is  too  large. 

This  part  of  the  computing  scheme  is  accurate  to  second 
order.   It  was  used  on  two  types  of  problems. 
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a)  Initial  value  problems  in  which  the  field  B(|,0)  (or 
pressure  N_^  +  e  N_ )  is  prescribed  at  T  =  0  and  all  other  data  kept 
constant  where  compatible   while  B(0,T)  is  kept  fixed. 

b)  Boundary  value  problems  (magnetic  piston)  in  which  all 
data  are  constant  initially  and  B(0,T)  is  prescribed. 

The  principal  object  of  a)  was  to  see  how  an  initial 
'discontinuity'  would  propagate.   If  our  equations  were  those  of 
magnetohydrodynamics  the  initial  discontinuity  would  be  resolved 
by  a  Riemann  shock  tube  configuration  with  a  shock  going  from  the 
higher  field  or  pressure  into  the  lower  field  or  pressure  region 
and  a  rarefaction  the  other  way. 

In  our  problems  this  change  is  prescribed  in  the  form 

B{0,i)    =  B    +  a   tanh  c{i  -  i    )  . 
'      o  o 

The  parameter  c  fixes  the  steepness  and  because  of  the  mesh  sizes 
involved  we  cannot  take  c  too  small.   The  changes  take  place  over 
approximately  half  the  Larmor  radius  of  the  ion.   For  zero 
pressure  the  field  first  flattens  and  then  generates  waves  which 
travel  ahead  as  is  expected  from  the  magnetic  piston  problem,  see 
Morton  [5].   If  the  amplitude  of  the  change  in  the  initial  B 
exceeds  about  1.4  the  wave  breaks  almost  immediately. 

This  pattern  is  repeated  if  pressure  is  present  but  the 
wave  breaks  for  far  less  increase  in  field. 


One  sees  from  (11)  and  (7)  we  can  either  prescribe  N  or  L 
constant  initially  but  not  both  if  B  ^  constant  initially. 
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In  problem  b)  a  change  in  strength  for  B  initially  constant 

is  prescribed  at  the  end  point.  For  compression  waves  the  end 

effects  were  very  strong  in  the  third  moments.   Here  we  describe 
only  rarefaction  results. 

6 .  Explanation  of  Corresponding  Figures 

Figures  1-5  show  how  the  wave  develops  in  the  B  field.   The 
dotted  line  is  the  initial  field.   The  solid  line  shows  B  at 
times  T  =  0,1,2,3,4.   Beyond  T  =  4  there  are  reflection  effects 
from  the  ends.   Figure  6  shows  how  the  ion  temperature,  initially 
.07,  develops  and  Figure  7  a  comparison  with  electron  pressure 
(adiabatic  formulas)  at  T  =  4.   The  ions  are  hotter  in  the  com- 
pression and  cooler  in  the  rarefaction.   All  later  figures  are  at 
T  =  4  to  show  maximum  effect.   Figure  8  shows  p~  -  1,  the  increase 
in  specific  volume.   Figure  9  shows  mean  velocity.   Note  the  sharp 
dip  in  this  case.   If  U  goes  to  zero  or  approximately  sound  speed 
the  wave  breaks.   Figure  10  shows  the  mean  velocity  in  the  y- 
direction.   Figure  11  the  off-diagonal  pressure  tensor  term  is  very 
small.   Figure  12  shows  how  little  the  pressure  in  the  y-direction 
changes  in  comparison  with  that  in  the  x-direction.   Figure  13 
shows  the  special  function  F.   Figure  l4  confirms  that  the  electric 
field  in  the  x-direction  is  0(l/e)  in  comparison  with  the  y-field 
in  Figure  15-   Note  there  is  a  considerable  end  effect  in  the 
Ep  =  E   field.   Figures  l6  and  17  show  the  third  moments  (see  the 

next  section  for  the  method).   The  largest  is  barely  of  the  order 

^/2 
(N  )  '^   and  shows  no  regular  pattern.   The  ends  for  a  length  l/2 
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have  been  smoothed  out.   The  condition  of  reflection  used  at  the 
ends  leads  to  very  bad  behavior  and  it  is  not  known  what  suitable 
relaxation  condition  should  be  made. 

Figures  l8-25  show  the  same  problem  with  half  the  rise  in 
the  field  but  N_^  'v  .1.   Qualitatively  the  same  phenomena  occur. 
However  the  third  moments  are  considerably  smaller  as  seen  in 
Figures  26  and  27,  certainly  negligible  for  any  self-consistent 
requirement.   Figure  28  shows  the  percentage  error  (maximum  2-4) 
in  the  density  computed  in  the  same  way  as  the  third  moments  and 
compared  with  the  solution  of  the  partial  difference  scheme. 
Figure  29  shows  that  in  one  place  the  particles  collected  and  in 
another  were  missing  —  there  were  80  in  each  mesh  box  to  start  with. 
The  two  points  correspond  with  the  two  minima  of  the  E  field,  see 
Figure  24. 

Figures  50-35  show  a  similar  Riemann  problem  with  an  initial 
step  in  pressure  in  a  constant  magnetic  field.   The  same  kind  of 
phenomena  occurs.   The  third  moments  however  develop  much  more 
shape.   Hopefully  this  problem  could  be  solved  self -consistently 
by  improving  the  mesh  size  by  a  factor  of  say  5- 

The  maximum  rarefaction  was  produced  by  cutting  B  to  half  its 
value.   Figures  36-4-5  illustrate  the  results  in  this  case.   Note 
that  the  pressures,  Figure  yj ,    may  be  taken  so  that  p  ~  1.   Note 

from  Figure  40  that  the  gas  density  has  been  halved  also.   The 

2 

third  moments  are  of  the  order  of  N  and  negligible.   One  notices 

that  N-,  p  has  relatively  little  noise  and  its  shape  may  be  valid. 
Figures  44  and  45  show  a  remarkable  uniformity  of  noise  in  the 
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density  and  the  peaks  in  the  displaced  particles  do  not  seem  to 
be  correlated  to  the  fields.  The  model  thus  appears  to  be  very 
effective  for  rarefaction  waves. 

The  stability  of  any  of  these  computations  has  not  been 
investigated.   For  the  purpose  of  computing  the  third  moments  such 
small  time  steps  were  taken  (k/h  =  .1)  that  no  apparent  difficul- 
ties arose  in  this  connection  and  certainly  h/k  was  large  compared 
to  all  characteristic  speeds. 

The  time  consuming  computation  is  involved  in  computing  the 
neglected  third  moments.   The  method  is  discussed  in  the  next 
section. 

7.  Computations  of  the  Third  Moments 

To  compute  the  third  moments  of  a  distribution  in  a  given 
field  we  follow  the  particles.   In  the  case  of  the  ions  this  means 
to  solve  the  transport  equation  (la)  by  finding  its  characteristics 
(the  particle  paths)  which  satisfy,  with  E  =  eE  , 

^  =  u 
dt 

^  =  E  +  evB 
dt 

5^  =  eE  -  euB  . 
dt     y 

This  was  done  by  the  Runge-Kutta  method  using  a  three  point  scheme. 

Since  E,  E  and  B  are  only  known  at  the  mesh  points  of  the 
Lagrangian  system  we  transformed  from  x,  t  to  ^,  T  and  interpolated 
the  fields  on  the  original  mesh. 
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Note  that  this  method  also  would  work  if  we  were  using  a 
completely  self -consistent  system  with  the  third  moments  present. 
Each  characteristic  is  computed  a  step  forward  after  the  field  at 
the  later  time  step  has  been  found. 

Let  us  assume  now  that  the  characteristics  can  be  computed 
with  an  arbitrary  degree  of  accuracy.   Then  (ll)  states  that  f 
is  constant  on  this  characteristic. 

At  T  =  0,  the  distribution  f   is  given  at  each  mesh  point 
and  at  each  mesh  point  in  some  region  of  u,v-space. 

To  each  mesh  point,  r,  s  in  the  initial  velocity  space  we 
assign  a  weight  factor  W   =  f (u(r, s ) , v(r, s ) )  and  an  area  factor 

I?  o 

a   which  together  represent  the  element  fdudv  corresponding  to 
that  mesh  point.   Then  initially,  if  the  density  is  1, 


(17)  y~  W   a   =  1  . 


The  initial  pressure  say  in  the  x-direction  is  represented  by 


(18)  y~  W   u^  a 


The  mesh  in  the  initial  velocity  space  was  taken  as  circular 
with,  for  example,  l6  divisions  of  the  angle  and  4  or  5  radii. 
The  largest  radius  difference  was  .1  but  depended  on  the  pressure. 
Increasing  the  number  of  points  to  l6o  in  the  basic  case  increased 
the  smoothness,  i.e.  the  fluctuations  were  cut  in  half.   This  is 
an  indication  that  in  the  future  with  better  codes  one  might  obtain 
smooth  enough  third  moments  for  self -consistent  computations  by 
increasing  the  number  of  mesh  points. 
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We  have  taken  f  either  as  a  step  function  or  as  a  Maxwellian 
in  which  case  the  density  p  =  1  normalized  the  Maxwellian  and  the 
exponent  was  determined  from  the  pressure. 

8.  Estimating  the  Third  Moments 

A  useful  lemma  for  estimating  is 

Lemma:   p~  fdudvd^  is  invariant  along  a  characteristic  in 
u, V, 4^  t-space . 

To  prove  this  one  integrates  the  differential  equation  (II) 
over  an  arbitrary  tube  jT  generated  by  a  set  of  characteristics. 
Thus,  applying  Green's  theorem  one  has 


/  f (t  +  ux  +  (E  +  evB)u  +e(E  -uB)v  )ds  =  0 
J^n    n  n    ^y      n 

where  t  ,  x  ,  u  ,  v  are  the  components  of  the  normal  on  the 
n'      n^   n   n 

boundary  hj    of  3".   The  factor  vanishes  on  the  characteristic 
boundary  and  choosing  h'J     initially  and  at  t  to  be  t  =  const,  we 
see  that  /  fdudvdx  is  independent  of  time.   But  /  fdudvdx 

r    -1  °' 

=  /  fp   dudvd|.   This  proves  the  lemma. 

a 

Next  we  note  that  any  moment 


G 


=  rf(u  -u)^(v  -  v)"^  dudv 


may  be  written  as 
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+H 

G  =  jjj    g(u,v,^,t)fdudvd^  +0(H^) 
-H 

+H 
=  17    pg(fr^)dudvd|  +  0(H^)  . 


We  can  evaluate  this  Integral  by  writing  it  as 
'  (pg)fp   dudvd^.   Here  a  is  a  region  of  u,v,|-space  con- 


n 

^n 
taining  a  particle  P  .   It  is  chosen   to  be  such  that  it  is  the 

end  of  the  characteristic  tube  at  t  generated  from  t  =  0  by  the 

three-dimensional  unit  o*   formed  by  the  mesh  unit  in  the  initial 

velocity  space  and  the  interval  ^*--^h,    i*  +  ^  h   where  i*   is  the 

initial  position  of  P  . 

n 

By  our  lemma. 


/  (pg)p   fdudvd^  =  /   pg(p   f)*du*dv*d^ 


^n 


=  r  pg(p~-^f  )*dr*d0''d|' 


* 

"n 


where  (p  ^f)*  is  the  value  of  p'-^f  at  r* ,    d* ,    C   with  r*,  0*  polar 

coordinates  in  the  initial  velocity  space.   This  integral  may  be 

-1   •)(■  "5 

approximated  by  (pg)^p  [i    )w^^ha   with  an  error  of  K^.   Here 

K  =  max  (hjK-,  )  where  K-,  is  the  maximum  mesh  width  in  r  ,9*-space, 

(pg)   Is  the  value  of  pg  at  the  point  P. 


Of  course  this  choice  cannot  be  made  for  particles  too  close  to 
the  boundary  |  =  H.   These  particles  will  be  denoted  by  Q^  and 
discussed  later. 
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The  particles  Q  would  produce  an  error  if  dropped  of  the 
order  h/H.   In  our  computations  where  h/H  =  l/4  we  followed  the 
following  method  which  smooths  out  the  end  effect  and  reduces  this 
error. 

The  method  depended  on  the  fact  that 

+H 
/fgdudv=/       fg    q(g')dudvd^'    ^Q(^2^ 

-H     y'q(?')d^' 

and  q  was  chosen  to  vanish  at  ?'  =±H  and  to  be  linear  on  (-H,0) 
and  (0,H).   Then  /  fg  dudv  was  approximated  by 

+H 

)         (Pgq)n  ^rs^rs^^  f       q^^')'^^'  ' 

-H 

(In  all  our  problems  p   (^  )  =  1.   Also  H  was  taken  to  be  2h. ) 

It  was  found  as  a  check  that  the  density  computed  this  way 
never  differed  from  that  computed  as  V~   in  the  p.d.e.  problem  by 
more  than  a  few  percent.   When  the  density  variation  was  as  much 
as  20  percent  the  maximum  error  was  6  percent.   The  third  moments 
were  at  their  largest  about  10"-^  which  was  the  same  order  as  K^. 
It  should  be  added  that  near  the  ends  of  the  total  ^  interval  the 
error  in  the  density  could  rise  as  high  as  20  percent  when  the 
density  had  hardly  changed. 

The  above  estimates  have  assumed  that  the  particle  paths  had 
negligible  error.   One  should  note  that  the  error  accumulates  in 
time  but  in  fact  this  will  be  of  lower  order  than  K^. 
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Appendix  I  -  Finite  Electron  Pressure 

For  finite  electron  pressure  the  speed  of  most  particles 
must  be  Od/e)  to  guarantee  that  N  =  e   /  f(u-u)  dudv  =  0(1)  and 
/  f  dudv  =  0(1).   We  make  the  change  of  scale  U  =  eu,  V  =  ev  and 
observe  that  the  transport  equation  for  the  electrons  becomes 

f  ,  +  Uf   -  (E+VB)f   +(eE  -UB)f   =  0 
-t     -X         '  _u   ^   y      -V 

where  we  have  set  E   =  E/e  in  accordance  with  (le). 

To  lowest  order  in  e  this  equation  reduces  to 

Uf   -  (E +VB)f   +  UBf   =  0  and  describes  the  transport  in  a  sheath 
-X   ^       -u     -v  ^ 

(i.e.   E  =0),  see  Sestero  [  2  ].   The  electrons  have  closed  paths 
in  (X,U,V) -space  which  satisfy 

U  f  =  _E  _VB 

dx 

X  X 

Setting  /   B  dx  =  A  we  have  then  with  ^  =  E  dx 

X*  -oo 

V  =  A+V""  ,    I^  +  (l)  +  4  +  V*A  =  ^+^* 


where  U*,  V  are  the  values  at  x  =  x*.   The  paths  are  closed 
since  A  -^ +oo  for  x  -»•  ±  oo  and  (|)  is  bounded.   The  particle  motion 


Note  that  the  full  equations  in  the  steady  state  require  every 
particle  to  go  from  -oo  to  +oo  ,  see  [9]' 


^7: 


given  by  the  Hamiltonian  has  an  ignorable  coordinate  to  lowest 
order  In  e.   Applying  the  theory  described  by  Gardner  [8]  we  have 
as  adlabatlc  Invariant  the  area  swept  out  by  an  orbit  In  the 
appropriate  phase  space.   This  area  a  Is  determined  by  Its  value 
at  X  =  X  .   The  distribution  f_  Is  therefore  to  lowest  order  a 
function  of  a  only.   This  problem,  while  theoretically  possible. 
Is  difficult  to  compute. 


Appendix  II 

Here  we  shall  show  that  for  a  distribution  of  ions  moving 
in  a  steady  field  we  have  the  second  moment  N  =  /  f(u-u)dudv 
behaving  like  p  ,  see  Chew,  Goldberger  and  Low  [6].   This  lends 
plausibility  to  the  calculations  with  neglected  third  moments  and 
small  off-diagonal  components  of  the  stress  tensor. 

Let  f  satisfy 

(1)      uf  +E(x)f   =  0   with   f  =  e~^^^"^^    for  x  =  ao  . 
x      ^  u 


Then 


2        /    2 
fu  dudv  -v  /  fu  du 


/    2 
But  the  element  udu  is  invariant  in  x.   Hence  /  fu  du  may  be  written 

as  /  f*u(x)u*du*  where  f*  is,  let  us  say,  the  value  of  f  at  co  for 

a  particle  with  u  =  u*  and  u(x)  is  its  velocity  at  x.   Then  since 

X 

^  =  u  J^  =  E  we  have  u^  =  2^   +  u*^  where  (f)  =  /   E  dx.   Thus 

oo 
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+00         _  p  1 

-GO 


where  we  assume  that  E  Is  such  that  there  are  no  closed  paths  in 
(x,y ) -space. 

We  set  k(u*-u)  =  |,  N  =  /  f(u-U)^du  where  U  =  /  fudu  //  fdu. 
By  the  conservation  of  mass  /  fudu  =  const,  which  we  normalize  to 
be  u  and  U  =  u  for  x  =  oo  .   Then 

2  ^ 

N+ uU  =  /  fu^du  =  ^  /  e~^  (  (u+^)  +  2<}))^(|  +u)d|  . 

We  expand  in  inverse  powers  of  k  and  find,  since  odd  terms  in  i 
vanish, 

N  +UU  =  k"^u(u^+  2<t))^  /  e"^  d|(l+|  -^ ^'k   — ^ 9  ^  +  .  ..) 

-'  2  u^+2(t)  k^   ^  (u^+  2(1))^  k^ 

1  -1  -1 

=  u(u^  +  2(t))^  +  :^  NQ(5U(u^  +  2<t))  ^  -u^(u^  +  2(t))  ^)  . 


But 


U  =.  u    f  r*{u*^+2^)~^   u^du* 


/      2  -^  ^"^ 

=  ku  (  f  e'^    ((u+|)^  +  2<i))  ^  (|  +  u)d? 

1  _  1  _2 

=  (u2+  2(1))^  +  I  N^(u2+  2(1))  ^  -  I  N^u^(u^  +  2(1))  ^ 


similarly.   Thus  we  find  with  pU  =  u  to  lowest  order 


if9 


_1 

N  ~  N  u^(u^+  2(j))  ^  +  0(^ 


i?' 


or 


N  ~  N  p^  +  O(N^) 


as  it  was  desired  to  show.   Note  that  the  neglected  term  is  of 
order  4  in  the  thermal  velocity. 
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Part  III  -  Shocks 


1.  The  Shock  Problem 

Observations  [J.  Paul  [I7]  etc.  [l8,19]]  of  transverse 
shocks  indicate  that  there  are  two  critical  shock  speeds  or  Mach 
numbers  when  the  observed  phenomena  for  transverse  shocks  change 
character.   Let  us  call  the  three  corresponding  zones  1,  2,  3: 

In  zone  1,  the  flow  is  oscillatory  with  wave  length  as  in  case  A, 

'    pe 
In  zone  2,  the  flow  has  a  sharp  initial  jump,  with  oscillations 

behind  of  wave  length,  c/oo 

°       pe 

In  zone  3?  the  flow  has  only  a  sharp  internal  jump,  i.e.  width  « 

c/co      .      Aside   from  this    jump    (  I       '  )       =  length  =  0    (c/od      ). 
pe  r>  pe 

This  description  is  fitted  by  Morton's  model  [3]  in  which  the 
three  zones  satisfy  the  following  criteria: 

Zone  1.  The  flow  is  subsonic  (using  gas  speed)  and  oscillatory. 

Zone  2.  The  flow  becomes  supersonic  in  an  internal  shock  but 
remains  oscillatory. 

Zone  3.  The  flow  has  no  oscillations  but  a  sharp  internal 
viscous  jump  to  supersonic  speeds. 

The  statements  describing  zone  2  could  be  altered  to  fit  the 
model  of  Morawetz  [6].   That  is,  by  introducing  enough  resistivity 
(nicknamed  anomalous)  the  internal  jump  can  be  dispensed  with  and 
replaced  by  a  strongly  damped  oscillation. 

Morton  and  Morawetz  are  based  on  a  gas  law  for  the  pressure 
tensors.   Modifications  of  such  models  to  treat  the  electrons 
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adlabatlcally  and  the  ions  by  the  electric  field  (see  case  A)  have 
been  made,  see  Marder  [7]  and  Part  II  of  this  paper.   The  basic 
observations  for  zone  1  remain  the  same. 

For  a  zone  with  a  jump  It  seems  reasonable  to  hypothesize  an 
Inner  layer  governed  either  by  large  charge  separation  fields  or 
by  turbulence  or  more  probably  by  both. 

2.  Turbulent  Inner  Shock 

In  a  purely  turbulent  model  of  an  Internal  shock  we  note  that 
even  If  the  flow  Is  highly  dependent  on  time  and  space  but  Is  con- 
fined to  a  narrow  region  any  equation  that  can  be  written  In  con- 
servation form  will  yield  a  possible  Jump  condition. 

We  examine  the  equation  for  U,  in  {J>1) ,    Part  I.   The  term 
pU*'VU*  may  be  written  as  dlv  pU*U*.   Thus  the  conservation  law  is 
to  be  found  by  Integrating  with  respect  to  time  and  space,  combining 
with  the  conservation  of  mass  and  obtaining 


pU  dx 


+    /         /    dlv    [pmj+M^+  e  Mg  +pU*U*]_dxdt 


/    (B-VB    -   ^igl     )    dxdt   =   0    . 


H 


The  space  integration  is  over  the  turbulent  layer  to  some  variable 
point.   Since  the  layer  is  supposed  thin  the  first  term  is  small  and 
the  only  contribution  from  integrating  dlv  terms  will  be  the  x- 
component.   Ordinary  one  fluid  theory  differs  from  this  only  in  the 
term  pU*U*  and  in  the  off  diagonal  terms  in  the  pressure  tensor. 
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Thus  in  an  average  sense  with  respect  to  time  and  y,z-space 
(indicated  by  bars),  the  average  momentum  is  conserved: 


B 


pU^U+(M^+  e  M2)^+pU^U*-  B^B+^^  k 


=  0 


where  k  is  a  unit  vector  in  the  x-direction  and  subscript  x  is  the 
x-component.   The  conservation  of  momentum  in  the  x-directlon  then 
,    yields  since  B  =  0 


P^x  +  ^^1  +  ^^^^2)x   +  pV^  +  ^/^ 


=  0  . 


To  obtain  an  effect  like  viscosity  there  would  have  to  be  a  term 
that  changed  sign  with  SU^/Sx.   This  leads  us  to  the  following  ansat2 

We  hypothesize  that  turbulence  is  produced  by  an  instability 
that  depends  on  the  magnitude  of  the  electric  field  in  the  x- 
direction  (i.e.  two  stream  instability  of  one  kind  or  another). 
Furthermore,  we  assume  that  the  average  effect  is  to  keep  B,  curl  B 
constant.   This  modifies  the  average  x-momentum  equations  to  the 
form 


pU^  + (M^  +  e^Mg )^  +  P  ^ ( curl  B )^  +  vE^ 


=  0 


where  [   ]  means  the  jump  from  say  x^  to  x,  and  (M^ +  e  M^)^  is  the 
scalar  pressure.   Here  v  is  some  proportionality  factor.   Since  the 
layer  is  assumed  thin  the  force  governing  all  motion  is  the 
electric  field. 
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Let  d(t)/dx  =  E  .   Then  all  flow  quantities,  p,  U  ,  M-,  ,  M„ 
can  be  seen  to  be  functions  of  ^.      The  momentum  equation  is  then, 
dropping  bars. 


V(})^  =  p(^) 


which  Integrates  to 


X  =  V 


f^  d(b 


J  p(()) 


where 


p  ( (f) )  =  pU^  +  (M^  +  ^^^2  ^x  +  P  ^  ^  ^^^1  ^  ^ 


2-M 

X 


with  M  the  total  momentum.   It  is  also  reasonable  to  assume  p((f))  >   0 
in  which  case  the  shock  layer  has  monotonically  increasing  cj)  in  the 
mean. 

Note  that  the  shock  width  is  proportional  to  v  and  roughly 
the  jump  in  potential  (|) .   The  state  ahead  is  given  by  a  zero  of 

p((t))  =  pU^  +(M^  +  e^M^)^  +p"^(curl  B)^  -M 

determined  by  one  upstream  state. 

One  hopefully  looks  for  a  second  root  to  find  the  state 
behind  the  shock.   For  this  a  necessary  condition  is  that  p'((l)) 
vanishes  somewhere.   Using  [6],  to  describe  the  ions  in  a  potential 
field  along  with  their  reflections  and  somewhat  crudely  treating 
the  electrons  for  this  argument  only  as  an  ordinary  gas,  p oC  p^ , 
one  finds,  since  by  [6], 


5^ 


^  (P^x+^lx)  =  4  /^^'^-^^   =  ^  /^^'  iu*^-2^du*dv* 


fu""  Ju*^  -  2([)  duMv*  =  -p  , 


that 


p'(i))  =  -p  +p^(c^  -p"^(curl  B)^) 

where  c_  Is  the  gas  sound  speed  of  the  electrons.   Thus  p'((j)) 

2 

will  vanish  if  c_  is  large  enough  and  pi  >•  0.   At  the  maximum 

value  of  (j),  again  by  [6],    pi  will  vanish  and  then  become  positive. 
If  it  becomes  large  enough  it  can  produce  a  second  root  to  p((}))  =  0, 
Thus  one  may  reasonably  anticipate  this  type  of  transition  for 
strong  enough  shocks. 

However,  independently  we  may  examine  the  states  behind  and 
in  front  of  this  special  layer  drawing  from  this  very  crude  model 
only  the  minimum  information.   We  proceed  to  do  this. 

3.  Transverse  Shock  Structure  for  Larger  Mach 
Numbers  with  an  Inner  Shock 

On  the  basis  of  the  model  described  in  [6]  and  the  layer 
described  in  §2,  we  hypothesize  certain  properties,  of  a  steady 
shock,  i.e.  at  rest  in  an  appropriate  frame.   The  flow,  which 
depends  on  x  only,  consists  of  four  regions: 

X  <  x  :   All  ions  with  u  negative  are  being  turned  around  by 
the  magnetic  field.   The  electric  field  is  0(1).   The  basic  width 
is  half  the  Larmor  radius  or  about  c/cd^. 
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x^   <  X  <  X-.:      The  ions  are  being  thrown  back  by  an  Increasing 
electric  field  0(l/e).   The  density  and  mean  x-velocity  are 
functions  of  ^,    and,  see  [6],  the  maximum  value  of  (j)  for  the  whole 
range  of  x.   Here  di/dx  =  eE_.   The  basic  length  is  c/co   ,  our 
basic  length.   The  electrons  are  essentially  adiabatlc.   All  we 
shall  need  to  use  is  that  the  mean  electron  velocities  are  given 
by  the  guiding  center. 

x-^   <  X  <   Xp:   A  black  box  in  which  we  assume  that  (j)  increases 
monotonically  and  across  which  B  and  curl  B  are  constant.   (it 
could  be  a  turbulent  boundary  layer  as  in  §2. ) 

Xp  **  X  <  GO  :   The  ions  are  again  governed  by  (j)  since  E=  0(l/e) 
The  density  is  a  function  of  (j)  determined  by  the  maximum  of  ^. 
Basic  length  c/oo   .   The  electrons  are  assumed  to  have  mean  veloci- 
ties  satisfying  the  appropriate  guiding  center  equations.   The 
electron  pressure  is  assumed  finite  here  so  this  is  a  bad  approxima- 
tion near  Xp.   However  for  large  x  where  the  fields  are  varying 
slowly  in  the  present  length  scale  it  is  a  valid  model  unless  the 
distribution  of  the  electrons  is  highly  skewed. 

Considerable  information  can  be  obtained  on  the  basis  of  this 
model,  in  particular,  the  final  values  of  pressure,  field  etc.  can 
be  determined  theoretically.   So  at  this  point  it  is  important  to 
restate  what  is  assumed;   (1)  The  nature  of  the  domain  x-,  <  x  <  Xp 
as  described  and  (2)  that  the  mean  velocities  of  the  electrons  in 
Xp  "^  X  *=  CO  satisfy  the  guiding  center  equations. 

Let  U,  V  be  the  (x,y)  velocities  of  the  guiding  center  motion. 
Then  from  II  say,  setting  E  =  E/e 
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^US   =   +Ey-UB,  0=|+VB 


Which  yields  eV  =  -E^/B.   From  curl  E  =  B^  =  0  we  have 
E  =  const.  =  1  in  our  present  normalization. 
Thus 

BU  =  1  -  U  d(E/B)/dx 


eV  =  E/B 


Reverting,  see  (3)  and  (4),  Part  I,  to  curl  B  =  J  instead  of 
using  (1),  Part  II,  we  find  ^  =  -£Vp_  =  -eVp  =  -em  ^  =  -  ^  since 
our  mass  flux  m  =  pU  is  constant  and  normalized  to  1.   Thus 

(1)  U  ^  =  I  ,    u  ;^  (|)  =  -1+UB 

dx   B  '      dx  ^B^ 


or 


d  rrr   dB. 


(2)  „^(trH|)  =  .1+UB 


It  is  easy  to  show  by  conservation  of  mass,  charge  and  charge 
neutrality  that  U  is  also  the  mean  x-velocity  of  the  ions  and  that 

U  =  p    is  a  function  of  (j).   Since  each  ion  with  x-velocity  u 

^..r,.         du     du   „   di  ,   -,     ^_         , 
satisfies  -r-r   =  u  -3—  =  E  =  -5-'-  to  lowest  order, 
dt     dx       dx 

The  actual  function  is  different  for  the  two  regions  under 

consideration,  i.e.   x  "<  x  <  x,  and  x„  <  x  <  00  .   We  denote  these 
'  o  12 

two  functions  by  R-,((j))  and  R^{^)    respectively.   It  is  technically 
impossible  to  determine  R-iCi")  explicitly.   The  governing  conditions 
are  described  in  Morawetz  [6].   For  our  purposes  it  suffices  to 
know  that  both  functions  are  uniquely  determined  by  giving  the 
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maximum  value  of  ^   provided  ^   is  monotonic  up  to  Xp.   Integrating 
(l),  from  Xp,  B  and  (f)  are  determined  as  functions  of  x  or 
B  =  B^{^,    \)^^,    B{x^))    for  X  >  Xg  or 

(3)  h   =  Mb)  . 


Note  also  that  B(x„)  =  B(x  ).   Thus  from  (2)  we  find,  with  U  ^  =  ii , 


dx 


dHp  =  (-l+UB)dB,  and  free  choice  of  a  constant. 


(4)  Jti^+H^(B)  =  0   for  X  >  Xp  . 


The  crux  of  our  argument  is  the  following.   If  the  solution 

behind  the  shock  is  to  be  eventually  constant,  it  must  tend  for 

dR 
large  x  to  satisfy  t]  =  -1  +  UB  =  0,  (i.e.   E  =  -U  x  B  and  -^^  =  0)  . 


dx 


In  other  words  it  must  pass  through 


dH 
Tl  =  0  ,    ^=   0  . 


In  the  neighborhood  of  this  point  in  the  B,  ri  plane 
-p  T]  +Hp(B)  =  0  must  be  a  pair  of  straight  lines.   (The  case  where 
the  curves  degenerate  to  a  point  corresponds  to  weak  shocks,  see 
[  3  ]  and  cannot  be  attained  without  resistance.) 

Thus  near  the  end  state  with  magnetic  field  B  =  B  we  must 
have 

(5)  H2(B)  <  0 

and  at  the  end  state 
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H2(B*)  =  0 
(6) 

-1+  u(b*)b''  =  0  . 


The  function  E^{B)    is  uniquely  determined  by  (4)  by  giving 

dR 
the  values  of  B,  i]   at  x  =  x^.   These  in  turn  are  given  by  B,  ^, 

dR 
and  U  at  X  =  Xp.   But  B,  -r—  are  constant  across  x-,  <  x  -^  Xp  and 

hence  are  determined  at  x,  by  (j)    .   U  _   Is  determined  by  R^{^) 

hence  by  (J5(xo)  and  <j)„„^  but  ({)(xp)  =  (f)^^^  provided  <{)  is  decreasing 

for  X  >•  x^.   Hence  B^  is  determined  by  cj)     from  the  equation 
(-■  max 

HpCB*)  =  0.   But  also  from  (6),  -1  +  U(B*)B*  =  0.   Hence  there  is 

one  condition  on  i    which  determines  it. 

^max 

Hence  the  final  value  of  B  is  determined  by  the  requirement 
that  we  end  in  a  constant  state.   (This  idea  appeared  first  in 
Morton  [3] . ) 

At  the  same  time  the  final  values  of  <[)  p,  U  etc.  are  all 
determined. 

We  examine  the  requirement  that  ^   is  decreasing  and  the 
meaning  of  the  condition  (5),  Hp(B)  <   0,    see  (4),  or  by  (l) 

B 

Hg  =  /   (-1+  UB)dB  =  -(B-  B"")  +  ^  -^*    . 
B* 

The  question  reduces  to: 


when  is 


B  -  (f) 


tB 

>•  0  ? 
JB* 


At  B  =  B*,  -T^  =  1  and  so  we  compute  — %.     If  this  quantity 
^^  dB^ 

is  negative  then  we  have  our  desired  type  of  endpoint. 
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From  the  momentum  equation,  (Part  II,  (l)  for  u  Independent 

2 

of  t),  pU  +  (M-,  +  e  Mp)   +  -p-  =  const,  we  can  determine  the  Jump 

across  the  "black  box"  in  electron  momentum.   We  write  the  total 

2 

pressure  M-,  +  e  Mp  as  p  and  regard  it,  as  we  may,  as  a  function  of 

X  and  hence  of  ^   and  hence  of  p.   We  denote 

dp  _   2 
dp   -    ^ 

Then  since  pU  =  1 

dU+  C^  dp  +BdB  =  0 

or 

du  -  C^  ^  +  BdB  =  0  . 
U 

2 

Hence  -^  =      ~    ^ .   But  UBdB  =  d^   by  the  current  equations.   Hence 
dB   y2  _  ^2 

dB^         ^^       U^-C"^ 

ah 


where  A  =  JB  /p  is  the  Alfven  speed  behind  the  shock.   Hence 

dB^ 
is  positive  if  the  flow  speed  is  less  than  gas  "sonic  speed"  as 

defined  here.   Lastly,  -^^^OinXp^x^oo  provided 

(C^-  U^-  A^)/(C^  -  U^)  is  positive.   This  will  hold  if  "the  black 

box"  represents  a  sufficiently  strong  internal  shock.   Equation  (4) 

guarantees  that  t5—  =  ?}  does  not  change  sign  so  that  ^  is  negative 

for  x„  <  X  <  oo  .   Hence  h  =   (t)(x„)  as  desired.   Thus  the  con- 

2  ^max    ^  ^  2 

figuration  agrees  with  the  conditions  used  to  describe  the  various 
zones.   The  chief  property  to  note  is  the  very  weak  assumptions 
about  the  motion.   The  worst  feature  is  of  course  the  neglect  of 
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the  gradient  of  the  electron  pressure  in  determining  the  current 
[cf.  (1)  and  the  equation  for  E,  equation  (l)  in  Part  II].  Note 
that  the  exponential  rate  of  approach  to  equilibrium  behind  the 


,   1    -Kx   .    .     V,   T^    /d  4) 
shock,  e    ,  IS  given  by  K  =  /  — 


dB 


B=B 


=  U/{A^+C^-U^)/(C^-U'^) 


B=B' 


Thus  the  width  of  the  final  region  scales  like  the  density  ratio 
for  very  large  Mach  numbers  since  U  ~  p~  . 
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This  report  was  prepared  as  an  account  of 
Government  sponsored  work.   Neither  the 
United  States,  nor  the  Commission,  nor  any 
person  acting  on  behalf  of  the  Commission: 

A.  Makes  any  warranty  or  representation, 
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in  this  report  may  not  infringe  privately 
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